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Results from the PHENIX experiment for the first RHIC run with Au-Au colhsions 
at ^/SNN = 130 GeV are presented. The systematic variation with centrahty of charged 
particle multiplicity, transverse energy, elliptic flow, identified particle spectra and yield 
ratios, and production of charged particles and tt^'s at high transverse momenta are 
presented. Results on two-pion correlations and electron spectra are also provided, along 
with a discussion of plans for the second run at RHIC. 



1. INTRODUCTION 

The PHENIX experiment has been designed to measure a broad variety of signals 
from both heavy ion and polarized proton-proton collisions at RHIC. The pursuit of 
penetrating probes generated in the early stages of the collision, combined with a program 
of hadron measurements, provides a detector with unparalleled capabilities to address 
observables sensitive to all stages of the collision process. This same detector is also very 
well-suited to the study of gluon and anti-quark contributions to the proton spinp|,||. 

The PHENIX detector consists of three spectrometers: two muon spectrometers cover- 
ing the full azimuth for 1.1 < \7]\ < 2.4 and a central spectrometer consisting of two arms 
each subtending 90° in azimuth and with \7]\ < 0.35. A central magnet provides an axial 
field, while each muon spectrometer contains a magnet that produces a roughly radial field. 
The central arms contain three tracking sub-systems: pad chambers (PC), drift cham- 
bers (DC) and time-expansion chambers (TEC); two forms of electromagnetic calorimetry 
(PbSc and PbGl); a time-of- flight hodoscope (TOP) and ring imaging Cerenkov counter 
(RICH). These sub-systems, together with a set of beam-beam counters (BBC) located in 
the region 3 < |?7| < 3.9, provide superb hadron and electron identification over a broad 
range of transverse momentum[Q. The muon spectrometers use cathode strip chambers 
in three stations for tracking (muTr), and five layers of larocci tubes interleaved with 
iron absorber for muon identification (muID). Global event characterization is achieved 
via a multiplicity and vertex detector (MVD) consisting of silicon strips and pads cov- 
ering |?7| < 2.5, and the RHIC-standard Zero-Degree Calorimeters (ZDCs), which detect 



4 



neutral particles emitted along the beam directions PJ^. The front end electronics for all 
sub-systems are clocked synchronously with the beam crossing frequency of 9.4 MHz. A 
set of Level-1 triggers derived from various sub-systems is used to initiate readout of the 
entire detector through a pipelined high bandwidth data acquisition system capable of 
archiving 20 MB/s. 



PHENIX Detector - First Year Physics Run PHENIX Detector - First Year Physics Run 




West Beam View East □ Installed □ Active 



Figure 1. Installed and active detectors for the RHIC Run-1 configuration of the PHENIX 
experiment. 

For the first physics run of RHIC in the summer of 2000, the portions of the PHENIX 
detector shown in Figure |I] were instrumented. Elements of all sub-systems, with the 
exception of the muTr, were in place and read out. Small subsets of the MVD and muID 
front end electronics were connected and tested as part of an engineering run. All other 
sub-systems were instrumented in fractions ranging from 25% to 100% of their ultimate 
aperture and were used in the physics results which are presented here. Independent 
minimum bias triggers were formed using coincidences between the BBC counters and 
between the ZDC counters. A total of approximately 5M events was recorded at ^/sNN = 
130 GeV. The primary trigger used for most of the results presented below is based on 
the BBC coincidence with an additional offline requirement that restricts the collision 
vertex to |z| < 20 cm. Whenever possible, physics quantities are presented as a function 
of centrality and/or the number of participants Npart- Details on the determination of 
these quantities are presented in the following section and in the associated references. 

2. GLOBAL OBSERVABLES 

The systematic variation of particle yields and the produced transverse energy with 
the number of participants reflects the underlying reaction mechanisms. For example, 
Gyulassy and WangjQ have emphasized that such a study can discriminate between cas- 
cade models and models which incorporate gluon saturation effects. PHENIX has studied 
the production near mid-rapidity of both charged particles and of transverse energy as a 
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function of centrality. The deposition of energy in the ZDC's is correlated with that of 
charge in the BBC's to provide an unambiguous mapping between these observables and 
the centrahty of the colhsion, as shown in Figure |. A Glauber model is then used to 
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Figure 2. BBC vs ZDC 
analog response, related to 
the centrahty classes of the 
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interval is the 0-5% cen- 
trality class, the next is 
5-10% and so on. 



determine the number of participants Np and the number of binary collisions Nc for each 
centrality class. (This methodology is used consistently for all such studies presented in 
this contribution.) 




Figure 3. a.) The pseudorapidity density per participant pair versus the number of 
participants, b.) The transverse energy pseudorapidity density per participant pair versus 
the number of participants. The band in each figure indicates the systematic errors. 



The number of produced charged particles is determined from the correlation between 
hits in two layers of Pad Chambers with the vertex location, which provides (on a sta- 
tistical basis) the charged particle multiplicity distribution in the interval \r]\ < 0.35. 
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After corrections for acceptance, efficiencies, decays and double hits, the charged particle 
pseudo-rapidity distribution dNch/drj is calculated for each centrality bin, scaled by the 
corresponding number of participant pairs, then plotted versus Np, as shown in Figure 
Also shown there are comparisons to the model predictions of EKRT P,p^ , which does 
not reproduce the trend of the data. (Results very similar to the PHENIX data were 
reported at this conference by the PHOBOS Collaboration [|TT].) 

A similar analysis has been performed for transverse energy measured in the PHENIX 



PbSc calorimeter |jT2| . A careful treatment of the contributions from produced energy 
in the aperture, the in-flux from scattering sources, and both the in-flux and out-flux 
from decays is performed to convert the deposited energy seen in the calorimeter to the 
equivalent transverse energy. The total transverse energy density per participant pair is 
then calculated in the same Np bins as used in the multiplicity analysis. The trend, shown 
in Figure is essentially identical to that found for charged multiplicity, and indicates 
that particle and transverse energy production is not simply proportional to the number 
of participants. The yield may be described as a superposition of terms proportional to 
participants and binary collisions. 

Yield = A-Np + B-Nc . (1) 

The results of this procedure are given in Table [l|. The extent to which the ratio B/A is 
significantly different from zero may be taken as evidence for the role of binary collisions 
in contributing to production of particles and transverse energy at RHIC[|13|. 



Table 1 



Quantity 


A 


B 


B/A 


dNch/di] 


0.88 ±0.28 


0.34^0.12 


0.38 ±0.19 


dEx / dr] 


0.80 ± 0.24 GeV 


0.23^0.09 GeV 


0.29 ±0.18 



The results of fits to Equation |I| for charged multiplicity and transverse energy near 
7] = 0. Note the dimensions of A and B are different for the two analyses. The =F symbol 
in the error for B is used to indicate that its error is largely anti-correlated with that of 
A. This anti-correlation is propagated in calculating the error in the ratio B/A. 



Also shown in Figure ^ are comparisons of both the charged multiplicity and the trans- 
verse energy per participant pair to the distributions measured by the WA98 collaboration 
at the CERN SPS|14], clearly demonstrating the increase of both quantities at RHIC. In 
the case of energy density calculated following the Bjorken prescription |]I5|, the value of 



4.6 GeV/fm^ is roughly 60% larger than that found at the SPS|T6|. It is interesting 
to note that this increase does not result from an increase in {dEx / drj) / {dNch/ drj) ; the 
value is very nearly equal to that found at the SPSfT^. 

Another global feature of hadronic production is the azimuthal pattern of emission, as 
parameterized by elliptic flow. Typically such analyses are performed with respect to the 
reaction plane determined for each event. It is also possible to determine the elliptic flow 
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Figure 4. V2 versus transverse momentum (a) and centrality (b). Errors are statistical 
only. 



pattern by measuring the auto-correlation function of particles in the event, that is, by 
calculating the correlation function C(A0), where A0 is the difference in azimuthal angle 
between two particles from the same event. PHENIX has performed such an analysis 
using charged particles with px > 200 MeV/c in the interval |?7| < 0.35. The correlation 
function is calculated using an event-mixing prescription to determine the background 
distribution, then fit to extract the (assumed positive) V2 coefficient: 

C(A0) oc 1 + 2|t;i|2cos(A0) + 2|t;2|^cos(2A0) (2) 

Standard PHENIX event characterization methods are used to study the dependence of 
V2 versus centrality in each bin. A sample of the results Ijl^ are shown in Figure ^. 
These trends, which are in good agreement with data from STAR[|^ and preliminary 
results from PHOBOS|Tl|, suggest that the high density matter formed at RHIC efficiently 



translates the initial spatial asymmetry into a corresponding one in momentum space. 



3. IDENTIFIED HADRONS 

The PHENIX central arm detectors have been designed to provide particle identification 
over the broadest possible momentum range |Q. The primary tool for charged hadron 
identification is the time-of-flight difference between the BBC and the highly-segmented 
TOF hodoscope, which spans A0 = 45" in the East spectrometer arm. The overall time 
resolution a ~ 115 ps permits unambiguous ir/K separation to at least pt = 1.5 GeV/c. 
Normalized minimum bias pr spectra for vr^, K^, p's and p's are shown in Figure ^ The 
shape of the spectra clearly depends on the particle species, with pions having the lowest 
(Pt) and protons and anti-protons having the largest. The dependence of the local slopes 
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Figure 5. Normalized minimum bias transverse momentum spectra for positive (left) and 
negative (right) identified particles. The error bars are composed of statistical errors 
and the systematic errors associated with acceptance and decay corrections. There is an 
additional 20% systematic error associated with the overall normalization. 



increases with both centrality and with particle mass, consistent with expectations from 



radial flow 19 



The ratio of anti-protons to protons at ?/ = is of particular interest, since it is a direct 
measure of the net baryon content in the central region. The p/p ratio has been studied as 
a function of both transverse momentum and centrality |2^, and found to be only weakly 
dependent on px and independent of centrality within systematic errors (Figure For 
minimum bias collisions, the p/p ratio in the interval 0.8 GeV/c < pr < 3.0 GeV/c is 
0.64 ± 0.01(stat.) ± 0.07(sys.). This value (which is consistent with that observed by the 
other RHIC experiments [p], pl| , p^ ) together with the measured yields shown in Figure ^ 
signifies that the central region in heavy ion collisions at RHIC is meson-dominated, as 
distinct from the baryon-dominated case at the SPS and lower energies, where the ratio 
is never greater than ~ 0.i p5| , pi . 

Charged pion pairs in the central arms have also been used to perform an HBT anal- 
ysis. Two separate analyses were performed, one using pions identified with the TOF as 
described above, and the second in which the 700 ps time-of-fiight resolution of the PbSc 
EmCal (EMC) is used to identify pions with transverse momentum pt < 0.7 GeV/c. In 



9 



a. 

I Q. 



0.8 



0.6 



0.4 



0.2 



> • 



i i 

^ i 



PHENIX preliminary 



0.5 1 1.5 2 2.5 3 3.5 

Pj [GeV/c] 



Q. 

I Q. 



0.8 



0.6 



0.4 



0.2 



n ' r 



PHENIX preliminary 



J I I I I \ , \ , I I \ , \ L 

b 50 100 150 200 250 300 350 400 

Npart 



Figure 6. The p/p ratio versus transverse momentum for minimum bias Au-Au colli- 
sions (left) and versus centrality (right). For both figures, the error bars represent that 
statistical error on the ratio; systematic errors are indicated by the horizontal bars. 



addition to the obvious utility of comparing two independent data sets for this analysis, 
the larger acceptance of the EMC substantially increases the number of available pairs 
(by a factor of 5). Correlation functions are calculated in the standard qrside, Qrout, QLong 
projections of the relative momentum for both it~^7t~^ and n^n^ pairs and are then fit over 
the full 3D phase space to the form 

C2{(lTside, QTOut, iLong) = 1 + ^ GW[~{lTsideRsideY — {iTOutRoutY ~ {.iLongRhongY] ■ (3) 

The results, presented in Figure ^ and Table ^ show little if any variation from values 
obtained at lower energies. Given the significantly higher multiplicities and densities 
observed at RHIC, this is somewhat puzzling, and perhaps indicates significant dynamic 
effects on the radii from strong transverse expansion driven by these higher densities. 
In this context, it should be noted that the values of the radii and A (for the TOF 
analysis), as well as the dependence on pair transverse momentumpSl are all consistent 



with measurements from the STAR collaboration! 29 . 



4. PARTICLE PRODUCTION AT HIGH TRANSVERSE MOMENTUM 

The high center-of-mass energy provided by RHIC, corresponding to values of ^/snn 
where hard scattering at the partonic level is observed in p-p and p-p collisions, offers 
the exciting possibility of using perturbative probes amenable to quantitative calculation 
to explore hot nuclear matter. A first step in this program is the measurement of the 
transverse momentum spectrum for charged particles and its variation with the number of 
participants in the collision. PHENIX has performed such an analysis using charged tracks 



10 



Au+Au - 2n 
PHENIX PRELIMINARY 




Figure 7. Correlation func- 
tion for TT^ pairs from the 
EMC analysis. The one- 
dimensional projections of the 
three-dimensional correlation 
function are averaged over the 
lowest 80 MeV/c in the other 
momentum differences. The 
fits are performed over the full 
WO 50 100 150 phase space; the curves in the 
q^^ (MeV/c) projections are averaged in the 
same manner as the data. 



Table 2 



Data Set 


Rxout (fm) 


Rxside (fm) 


Rlong (fm) 


A 


EMC TT+TT+ 


4.4 ±0.2 


5.1 ±0.6 


5.9 ±0.4 


0.27 ± .02 


TOP 7r+7r+ 


6.2 ±0.5 


7.9 ± 1.1 


4.0 ± 1.2 


0.49 ± .05 


EMC TT-TT- 


5.1 ±0.2 


5.0 ±0.6 


5.9 ±0.4 


0.30 ± .02 


TOP n-TT- 


5.5 ±0.5 


5.8 ± 1.5 


6.7 ±0.9 


0.49 ± .06 



Results of the Bertsch-Pratt fits to the identical pion pairs in the EMC and TOP analyses. 
Only statistical errors are shown; current systematic uncertainties are < 1 fm. 



reconstructed with the drift and pad chambers of the central arm||26||. The raw distribution 
of tracks is formed for six exclusive centrality classes, again using the standard PHENIX 
event classification described in Section Corrections are made in each centrality bin for 
the spectral distortions due to momentum resolution and to backgrounds from scattering 
and decays. Results are shown in Pigure ^ for the interval 0.4 < < 5 GeV/c, which 
spans six decades in yield. 

In nucleon-nucleon collisions, production of high transverse momentum particles results 
from hard scattering at the partonic level. Since these cross sections are small, their 
contribution in nucleus-nucleus collisions, absent collective effects, should scale as the 
total number of binary collisions. Por the most peripheral (80-92% in centrality) bin 
shown in Pigure the shape and yield above pt > 2 GeV/c are in good agreement 
with an interpolation of hadron-hadron data|^ scaled by the number of binary collisions 
{{Ncoii) = 3.7 ± 2) calculated for this centrality bin. This is not the case for the most 
central collisions, which show a deficit at large pt with respect to Ncou scaling. This is 
presented in Pigure |^, where the scaled ratio 

YieldiCentral) / {NcouiCentral)) 



Yield{Peripheral) / {Ncoii{Peripheral)) 
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Figure 8. Normalized invariant 
transverse momentum spectra in 
the interval |?7| < 0.35, for six cen- 
trality classes. The curve compares 
the most peripheral class (80-92%) 
to the yield expected from p-p col- 
lisions scaled by the corresponding 
number of binary collisions. 



Figure 9. The pt dependence of the ra- 
tio of yields for central (0-5%) collisions to 
yields for two peripheral classes (60-80% and 
80-92%), scaled by the corresponding mean 
number of binary collisions for each central- 
ity class. 



is plotted for two different values of the peripheral reference set. The suppression at 
Pt > 3 GeV/c is inconsistent both with enhancements expected in that region for the 
Cronin effect ||28| and also seen in PHENIX data for mid-peripheral data when scaled 
to the p-p reference distribution!^. Similar results for such a suppression pattern in 
unidentified charged particles have been reported by the STAR collaboration |2^ . 



An independent analysis by PHENIX provides further insight into the nature of the 
suppression ||30||. The superb segmentation Ar] x = 0.01 x 0.01 and excellent resolution 



8.2%/ \J E (GeV) © 1.9% of the PbSc electromagnetic calorimeter are used to extract the 
vr" Pt spectrum by reconstruction (on a statistical basis) of their principal vr'^ 77 decay 
mode. Extensive studies based on the mixing of single photon showers into both real and 
simulated events are used to determine the variation of the reconstruction efficiency with 
multiplicity and of the background with the event multiplicity. A detailed simulation 
is used to estimate the contributions from background and decay particles. Results are 
shown in Figure |l^ for both a peripheral (60-80%) and central (0-10%) sample. As in 
the unidentified charged particle analysis, the peripheral data are well described with an 
interpolation of existing hadron-hadron data|^ scaled by the number of binary collisions. 
(An additional factor derived from p-p data of 7r°/(/i+ -|- h~) = 1/3.2 is used to scale the 
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Peripheral (60-80%) and central (10%) Jt° data 
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Figure 10. Normalized px spectra for 
7r°'s from peripheral (60-80%) and central 
(0-10%) Au-Au collisions. The bands com- 
pare the data to the yield expected from p-p 
collisions scaled by the corresponding num- 
ber of binary collisions; their width indicates 
the systematic uncertainty in the number of 
collisions used for scaling. 



Figure 11. The pt dependence of the ratio 
of 7r° yields for central (0-10%) collisions 
to yields for a peripheral collisions (60- 
80%), scaled by the corresponding mean 
number of binary collisions for each cen- 
trality class. The solid line represents the 
systematic error on the yields; the dashed 
line the additional systematic error from 
the uncertainty in the number of collisions 
used for scaling. 



measured yield of unidentified charged particles to neutral pions.) Once again, the central 
data fall well below the corresponding yield expected if high pt particle production scaled 
with the number of binary collisions. This is strikingly illustrated by the scaled ratio of 
Equation ^ for the central to peripheral yields, as shown in Figure |Tl]. 

A comparison of the ratio in Figure ^ for unidentified charged particles to that of 
Figure [Tl| for vr^'s suggests that the suppression may be more pronounced in the case of 
identified (neutral) pions. Should this be established (the present understanding of the 
systematic errors does not permit a definite conclusion) a consistent description would 
imply that charged pions are also preferentially suppressed as compared to the unidentified 
charged particles. Hints of precisely this behavior are seen in Figure |], where the trends 
in identified particle spectra suggest that the contributions from protons and anti-protons 
become comparable to those from pions for pt > 2 GeV/c. Further work is required to 
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improve both statistical and especially systematic errors to determine if the agreement 
of these data[^ with "jet quenching" predictions ||31| is in fact evidence for an enhanced 
energy loss mechanism in hot nuclear matter. Clearly essential to that program is a 
set of detailed measurements of both p-p reference data and proton-nucleus collisions to 
determine the quantitative value of the Cronin effect at RHIC energies. 



5. FUTURE PLANS 



A major component of the PHENIX 
physics program is dedicated to measure- 
ment of leptonic signals. For example, 
the 7r° spectra described in Section ^, 
while intriguing and valuable in them- 
selves, are also a prelude to the mea- 
surement of direct photons. In the cen- 
tral arms, virtual photons and vector 
mesons can be detected via their decays 
to e^e~ pairs, while open charm and bot- 
tom production are expected to domi- 
nate the production of single electrons for 
Pt >2 GeV/c. Measurement of these sig- 
nals requires n/e rejection in excess of 10'^ 
and very careful control of background 
contributions. Figure |I2| shows the first 
efforts in this program|^2|. Charged par- 



ticles are tracked using the PHENIX drift 
and pad chambers. Electrons are selected 
by requiring that at least three PMT's fire 
in the RICH, and are further identified 
via a tight matching cut between track- 
ing momentum and energy deposition in 
the EmCal. 
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Figure 12. The inclusive spectrum 
of unidentified hadrons, vr^, vr^'s 
and measured in Au-Au mini- 
mum bias collisions. 



The power of this combined particle identification from the various PHENIX sub- 
systems is apparent from Figure |T2|, where a clean electron spectrum is extracted at a level 
2-3 orders of magnitude below that of all charged particles. This approach will form the 
basis for future PHENIX measurements of vector mesons and open charm at RHIC, and 
is also applicable to photon measurements via external conversion. The measurement of 
electrons along with all other PHENIX analyses will benefit greatly from the significantly 
enhancements made to the detector for Run-2 at RHIC. As shown in Figure |13|, the com- 
plete aperture of the central arms will be available (compare to Figure |lD, much of the 
MVD will be instrumented, and an entirely new spectrometer to measure muons will be 
deployed. These additions in aperture and in capability, coupled with significant upgrades 
to the data acquisition and triggering system, should result in a hundredfold or more in- 
crease in the event sample obtained from Run-2. This will allow PHENIX to explore 
new signals and observables, to greatly increase statistical precision and understanding 
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PHENIX Detector - Second Year Physics Run PHENIX Detector - Second Year Physics Run 




West Beam View East □ Installed □ Active 



Figure 13. Installed and active detectors for the RHIC Run-2 configuration of the 
PHENIX experiment. 



of systematic errors on existing analyses, to obtain the vital p-p comparison data, and to 
begin a program of measurements with polarized protons dedicated to understanding the 
proton spin. 
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